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ABSTRACT 

Aims. We investigate the structure of the circumstellar disk of the T Tauri star S CrA N and test whether the observations agree with 
the standard picture proposed for Herbig Ae stars. 
1*%/ Methods. Our observations were carried out with the VLTI/AMBER instrument in the H and K bands with the low spectral resolution 

mode. For the interpretation of our near-infrared AMBER and archival mid-infrared MIDI visibilities, we employed both geometric 
and temperature-gradient models. 

^—1 Results. To characterize the disk size, we first fitted geometric models consisting of a stellar point source, a ring-shaped disk, and a 

Clj halo structure to the visibilities. In the H and K bands, we measured ring-fit radii of 0.73+0.03 mas (corresponding to 0.095±0.018 AU 

for a distance of 130 pc) and 0.85 ± 0.07 mas (0.1 1 1 + 0.026 AU), respectively. This /f-band radius is approximately two times larger 
than the dust sublimation radius of a;0.05 AU expected for a dust sublimation temperature of 1500 K and gray dust opacities, but 
approximately agrees with the prediction of models including backwarming (namely a radius of r;0.12 AU). The derived temperature- 
gradient models suggest that the disk is approximately face-on consisting of two disk components with a gap between star and disk. 
The inner disk component has a temperature close to the dust sublimation temperature and a quite narrow intensity distribution with 
a radial extension from 0. 1 1 AU to 0. 14 AU. 

Conclusions. Both our geometric and temperature-gradient models suggest that the T Tauri star S CrA N is surrounded by a circum- 
stellar disk that is truncated at an inner radius of «0. 1 1 AU. The narrow extension of the inner temperature-gradient disk component 
| implies that there is a hot inner rim. 

^2 ' Key words. Stars: individual: S CrA N - Stars: pre-main sequence - Stars: circumstellar matter - Protoplanetary disks - Accretion, 

, accretion disks - Techniques: interferometric 

l> 

■ 1 . Introduction surrounded by an additional extended halo of scattered starlight, 

which influences the precise determination of the disk size 
Near- and mid-infrared interferometry is able to probe the in- (Pinte et al 2008) The positions of observed TTS in the size- 
ner regions of the circumstellar disks of young stellar objects luminosity re i ation (Eisner et al. 2007) suggest, that TTS have 
(YSO) with unprecedented spatial resolution. However, the de- slightly larger inner disk radii than expected . However, if one 
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tailed structure of the inner gas and dust disks is not yet well- compares the TTS radii with predictions of models including 

g . known. In particular, the disks of T Tauri stars (TTS) are difficult backwarming (Millan-Gabet et al. 2007; Dullemond & Monnier 

■ - " to study because of their lower apparent brightnesses and the dif- 2010) the discrepancy disappears 

ficulty in spatially resolving them (e.g. Akeson et al. 2005a,b). , . . . , . „ ,. . „ , 

It is not yet known, for example, whether there is a puffed- *V^ S FT*' ,* e u mve f § ate the circumstellar disk of the 

up inner rim (PUIR) at the inner edge of TTS, as observed in ™ A S 9"\ N ' whlch 18 the more masslve , star m „ th ? J? 1 "^ 

several Herbig Ae/Be disks (Natta et al. 2001; Dullemond et al. f D CrA ' binary separation is approximately 1.4 (.150 AU) 

2001; Muzerolle et al. 2003; Monnier et al. 2005; Cieza et al. (Re ! pU ^ & - Zl " ne o ck ^1993; Ghez et al 1997) and its position 

2005). Furthermore, observations suggest that several TTS are an 8 le < PA > 18 *f < Ghez et aL WT). ™ e blnar y c ™P™f 

are coeval and have an age of «3 Myr (Prato et al. 2003). The 

properties of both stars are listed in Table 1 . 



* Based on observations made with ESO telescopes at the La 

Silla Paranal Observatory under program IDs 081.C-0272(A), 083.C- S CrA N ls a c l ass i cal TT S (McCabe et al. 2006). Its m- 

0236(C) frared excess suggests the presence of a dusty disk. The precise 

** Member of the International Max Planck Research School (IMPRS) determination of its spectral type is difficult owing to a strong 

for Astronomy and Astrophysics at the Universities of Bonn and veiling of the absorption lines (Bonsack 1961). McCabe et al. 

Cologne (2006) inferred a spectral type of K3 and Carmona et al. (2007) 
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Table 1. Properties of the S CrA binary components 
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Parameter 


S CrA N 


SCrAS 


60 


spectral type 


K3 


M0 


AO 


M» [Ma] 


1.5 ±0.2 


0.6 ± 0.2 


T, [K] 


4800 ± 400 


3800 ± 400 


20 


L, [L s ] 


2.30 ± 0.70 


0.76 ± 0.24 


Bry [10~ 16 Wm- 2 ] 


4.23 ±1.10 


1.59 ±0.57 


* 


mj [mag] 


8.6 


9.4 


> 


m H [mag] 


7.5 


8.3 


20 


m K [mag] 


6.6 


7.3 


distance [pc] 


130 


±20 


40 


A v 


2. 


¥' 




binary sep. ["] 


1.30±0.05'71-4 C 


60 



Notes. If not mentioned otherwise, the values are taken from Prato et al. 
(2003). For the spectral type of the primary, Herbig & Bell (1988) found 
K6 and Carmona et al. (2007) G5Ve. Ortiz et al. (2010) found a distance 
of 138± 16 pc, which is consistent with the table value and the one exten- 
sively discussed by Neuhauser & Forbrich (2008) (130 pc). 
Other references: <fl) Patten (1998), {b} McCabe et al. (2006), (c) Ghez et al. 
(1997) 
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AMBER 2008-06-09 
AMBER 2009-05-21 
MIDI 2005-05-30 
MIDI 2005-06-27 



Fig. 1. The uv coverage of the AMBER (red and green dots) and 
MIDI (blue and pink squares) observations used for modeling. 



obtained a spectral type of G5Ve. The veiling as well as the de- 
tection of a strong Bry flux suggest the presence of an accretion 
disk (Prato et al. 2003; McCabe etal. 2006). Schegerer et al. 
(2009) resolved the disk of S CrA N in the mid-infrared with 
VLTI/MIDI and modeled the spectral energy distribution (SED) 
and visibilities with the Monte Carlo code MC3D to constrain 
several disk parameters. 

The S CrA system is probably connected with Herbig-Haro 
objects: HH 82A and B are oriented towards a position angle of 
*95°, whereas HH 729 A, B and C lie in the direction of *115° 
(Reipurth & Graham 1988). The different PA can, for example, 
be explained by either two independent outflows from each of 
the binary components of S CrA or by regarding the HH ob- 
jects as the edges of an outflow cavity (Wang et al. 2004). The 
determination of the orientation of the circumstellar disk might 
clarify the exact relationship. 

Walter & Miner (2005) found that the secondary, S CrA S, 
can be brighter in the optical than the primary for up to one third 
of the time and that S CrA has one of the most rapidly varying 
brightnesses of the TTS. This variability is discussed in Graham 
(1992), who proposes that it is caused by geometrical obscura- 
tion as well as accretion processes and emphasizes the probabil- 
ity of clumpy accretion. 

Only very recently, Ortiz et al. (2010) observed light echoes 
in the reflection nebula around S CrA and suggested that 
the structure is similar to the Oort cloud in our solar sys- 
tem. With their method, they also determined the distance of 
S CrA (138+16 pc) confirming the previously derived distance 
of 130 pc (Prato et al. 2003; Neuhauser & Forbrich 2008). 

In this paper, we present VLTI/AMBER observations and 
the temperature-gradient modeling of the circumstellar disk of 
S CrA N. In Sect. 2, we describe the observations and data re- 
duction. In Sect. 3, we describe our modeling, which is discussed 
in Sect. 4. 

2. Observation and data reduction 

Our observations (program IDs 081.C-0272(A) and 083. C- 
0236(C)) were carried out with the near-infrared interferome- 
try instrument AMBER (Petrov et al. 2007) of the Very Large 
Telescope Interferometer (VLTI) in the low-spectral-resolution 
mode (R = 30). The observational parameters of our data sets 



are listed in Table 2 and the uv coverage is shown in Fig. 1. 
All data were recorded without using the FINITO fringe tracker. 
The observations were performed with projected baselines in the 
range from 16 m to 72 m. 

For data reduction, we used amdlib-3.0 1 . A fraction of the 
interferograms were of low quality. We therefore selected the 
20% data with the highest fringe signal-to-noise-ratio of both the 
target and the calibrator interferograms to improve the visibility 
calibration (Tatulli et al. 2007). Furthermore, it was impossible 
to reduce the H band data of data set I. We applied an equalisa- 
tion of the optical path difference histograms of calibrator and 
target to improve the visibility calibration (Kreplin et al. 2012). 

The derived visibilities are shown in Figs. 2 and 4. The mea- 
sured closure phase (Fig. 3) is approximately zero, which sug- 
gests that it is an approximately symmetric object. However, 
small closure phases are also expected because the object is only 
partially resolved. 



3. Modeling 

3. 1 . Geometric models 

To measure the characteristic size of the circumstellar environ- 
ment in the H- and /T-bands, we fitted geometric models to the 
visibilities. Figure 2 shows the observed visibilities and the fit- 
ted models, which consist of a ring (the ring width being 20% 
of the inner radius; Eisner et al. 2003; Monnier et al. 2005), a 
stellar point source, and an extended fully resolved halo (in only 
the lower two panels). This extended halo is assumed to rep- 
resent the stellar light scattered off a large-scale circumstellar 
structure (Akeson et al. 2005b). The motivation for a ring fit 
is that the model intensity distribution is expected to have a 
dominant PUIR brightness (Natta et al. 2001; Dullemond et al. 
2001). Ring-fit radii are often used in the literature to char- 
acterize the disk size of YSO and to discuss their location in 
the size-luminosity relation (e.g. Monnier & Millan-Gabet 2002; 
Dullemond & Monnier 2010). 

The flux contribution / star + /halo from the star (/ sta r) plus a 
halo of scattered starlight (/h a i ) was derived from the SED fit in 
Fig. 4 and amounts to 0.22 of the total flux (i.e. /star + /halo + /disk) 



1 http : //www . jmmc . fr/ data_processing_amber . htm 
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Table 2. AMBER observations 



Data set 


Night 


Bproj [m] 


PA[°] 


Seeing [" 


Airmass 


DIT [ms] 


Calibrator 


Calibrator diameter [mas] 


I 


2008-06-09 


16/32/47 


240/240/240 


0.91 


1.07 


100 


HD 183 925 


1.44+0.02" 


n 


2009-05-21 


60/69/72 


145/263/196 


0.98 


1.06 


150 


HD 170773 


0.38±0.026 " 



Notes. The parameters of the MIDI observations used for the modeling in Sect. 3.2 are described in Schegerer et al. (2009) and the uv coverage is 
shown in Fig. 1. References: (fl) Richichi et al. (2005), M Kharchenko & Roeser (2009) 




Fig. 2. Top: //-band (left) and Zf-band (right) visibilities of S CrA N together with ring-fit models (described in Sect. 3.1) consisting 
of a circular symmetric ring (ring width = 20% of inner ring radius r r ; ng; i n ) and an unresolved stellar source. Bottom: Models 
consisting of the same star+ring model as above, plus a fully resolved halo component. 



in the Zf band (2.2 pm) and 0.40 in the //band (1.6 pm). The total 
visibility V can be described by 



W = 1(1 _ /star - /halo)Vdisk + /starKtar + /halo Vha 



(1) 



where V stal = 1 and y ha io = 0. 

In fitting the data, we averaged the visibility data in each of 
the two bands. The H band visibilities were averaged over the 
wavelength range of ss 1.55-1. 75 pm and the Zf band visibilities 
over «1. 95-2.40 pm. 

For the inner fit radius rring.in, we obtained 0.73 + 0.03 mas 
in the //band and 0.85 + 0.07 mas in the Zf band. The linear radii 
in AU for a distance of 130 pc (see Table 1) are listed in Table 3. 

To test whether an elongation of the object can be measured, 
we also fitted an elliptic ring model. However, we were unable 
to detect any significant elongation. Because of the errors in the 
data, we can only derive a rough estimate of the elongation of 
< 30%. Therefore, we assume in the following modeling that 
the disk is oriented approximately face-on. 

3.2. Temperature-gradient disk model 

In this section, we present temperature-gradient models to help 
us interpret our observations. Each temperature-gradient model 
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Fig. 3. Closure phases measured with AMBER in the H and K 
bands. 



is the sum of many rings that emit blackbody radiation with tem- 
peratures T(r). For the temperature distribution, a power law is 
assumed (Hillenbrand et al. 1992), of T(r) = T (r/r T q . Here, 
Tq is the effective temperature at a reference radius ro. The pa- 
rameter q depends on the disk morphology and is predicted to be 
0.50 for flared irradiated disks and either 0.75 for standard vis- 
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Table 3. Ring-fit parameters in the H and K bands. The error in the distance measurement (130 + 20 pc, see Table 1) is included in 
the errors in the radii. 



Model 


Band 


' ring. in 


/halo 


/. 




/disk 






[AU] 










ring-star 


H 


0.120 + 0.020 





0.40 + 


0.09 


0.60 ± 0.09 


ring-star 


K 


0.134 + 0.027 





0.22 + 


0.05 


0.78 ± 0.05 


ring-star-halo 


H 


0.095 ± 0.018 


0.08 ± 0.05 


0.32 + 


0.14 


0.60 ± 0.09 


ring-star-halo 


K 


0.111+0.026 


0.06 ± 0.04 


0.16 + 


0.09 


0.78 ± 0.05 



Table 4. Overview of the parameter space scanned for the temperature-gradient models and parameters (with errors) of the best-fit 
models for one temperature-gradient disk (A, B) and two temperature-gradient disk components (C, D, E). For details about the 
individual models, we refer the reader to Sect. 3.2. 



Parameter 


Scan range 


A 


B 


C 


D 


E 


Constraints 








'out.l = r m,2 

T ut,i — 7i„2 


r out,l = r in,2 




r mA [AU] 


0.01-0.3 






11 +U06 

11 -0.04 


14 +0U3 


13 +UU4 

"' iJ -0.05 


fout,l [AU] 


0.02^1 






14+ 005 

"•"-0.03 


17 +009 

1 ' -0.05 


17+ 009 

1 ' -0.05 


n„.2 [AU] 


0.01-0.6 


0.05 + 0.01 


0.05 ± 0.03 


14 +005 

u - -0.03 


17 +009 


14 +004 

iH -0.05 


r om ,2 [AU] 


1-50 


> 4 


> 5 


> 26 


> 26 


> 26 


ZVi [T] 


100-3000 






1690+ 170 


1530!|» 


1500+^o 


T m2 [T] 


200-3000 


1900+ 90 


1890+?° 


660+ 90 

UU -160 


633+ 130 

UJ -120 


660+ 90 


?i 


0.35-0.85 






0.2 ± 0.2 


0.5 ± 0.4 


0.5 ± 0.4 


92 


0.3-0.9 


75 +0 1 


0.75 + 0.1 


0.5 + 0.1 


5 +012 

"• J -0.04 


5 +012 

""'-0.04 


/halo 


0-1 




0.2 + 0.1 








2 




16.9 


11.5 


3.2 


3.0 


3.0 



cous disks or flat irradiated disks (Chiang & Goldreich 1997). 
Wavelength-dependent visibilities and fluxes of these model 
disks depend on the disk inclination, the inner radius ri n , the 
outer radius r out , the temperature T m at r m , and the temperature 
power-law index q. As explained in Sect. 3.1, we assumed that 
the disk is viewed face-on. 

To determine the best temperature-gradient model, we calcu- 
lated the model disks (one- and two-component structures) for 
all combinations of the parameters with scan ranges described 
in Table 4 (s;400000 models). The total is a sum over the 
X^ ed of all visibility points for all six AMBER baselines and two 
MIDI baselines (see Schegerer et al. 2009), and the SED. We 
computed all combinations by running the model for six to ten 
steps per parameter and then chose the areas with the smallest 
value of x^ ed to obtain a finer mesh. We took into account only 
wavelengths with A < 20fim. The errors in the best-fit model 
parameters are 1-cr errors. 

The model assumes a stellar point source with the parameters 
T t = 4800 K, L, = 2.3L , distance = 130 pc, and A v = 2.8 
(see references for the stellar parameters in Table 1). The model 
parameters of all derived models (disk plus star, disk plus star 
plus halo, as well as several two-component disks) are listed in 
Table 4. 

3.2.1. One-component temperature-gradient disk models 

The simplest temperature-gradient model A (see Table 4) con- 
sists of only one single disk component plus the star. It cannot 
reproduce the SED as well as the MIR and NIR visibilities si- 
multaneously (xLa ~17). 

For the star-disk-halo model B, we added a fully resolved 
halo (V = 0) to the star-disk model A. This halo intensity 
distribution is assumed to represent the stellar light scattered 



off the large-scale circumstellar material and to have approxi- 
mately the same SED as the star itself (cf. Akeson et al. 2005b). 
Nevertheless, the improvement achieved was very small since 
the^ ed became only 11.5 (see Table 4). We also computed 20 
star-disk-halo models with resolved halos of different widths in 
the range from 0.05 AU to 100 AU, but the^ ed did not signifi- 
cantly improve since it remained larger than sj11.5 in all cases. 
As in the case of model A without a halo, it was only possible 
to fit either the SED and the AMBER measurements or the SED 
and the MIDI measurements. This suggests that a more compli- 
cated disk structure should be considered. Therefore, we intro- 
duced a second disk component in the following, but - as we 
preferred to adopt a two-disk modeling with a minimum number 
of free parameters - we omitted the weak halo. 

3.2.2. Two-component temperature-gradient disk models 

Model C is the result of our attempt to find a two-disk compo- 
nent model with the least number of parameters (Table 4). We 
therefore introduced constraints on the inner and outer radii and 
temperatures of disks 1 and 2, of T out j = T m ,2 and r out j = r m ,i- 
We obtained yfc. = 3.2 (the power-law index q\ is determined 
from the temperature slope between and r outj i and is there- 
fore no longer a free parameter). 

We also tested models with more parameters: In model D, 
'"out.i = r-m 2 is the only constraint and^ ed = 3.0. For model E, 
there are no constraints, but we still get a similar (3.0). 
Table 4 shows that only models C to E have &xLa fr° m 3.0 to 
3.2 and among those, model C has the advantage that it has the 
smallest number of parameters. 

Fig. 4 shows the intensity distribution of model C, together 
with all S CrA N observations. An extended ring is located be- 
tween 0.14 AU and 26 AU, whereas the inner ring extends radi- 
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Fig. 4. Temperature-gradient model C of S CrA N - Two top rows: NIR visibility. The model that fits visibility and SED simulta- 
neously (model C in Table 4) is indicated with the solid red lines. The six figures show the modeling of the data sets I (top) and 
II (bottom), described in Table 2. Third row: MIDI visibility. The two figures show the comparison between our model (red line) 
and MIR data adopted from Schegerer et al. (2009). Bottom left: Spectral energy distribution. To model the SED of S CrA N, we 
collected values from the literature and binned them (black dots are the dereddened SED points, see Appendix A). The green line 
represents the MIDI spectrum of Schegerer et al. (2009). The SED of model C is indicated with the blue line. It consists of the 
stellar contribution (Kurucz model, black dashed line) plus the contribution of two additional ring-like structures (black solid and 
dash-dotted line). Bottom right: Intensity distribution of the temperature-gradient disk of model C. The dashed white ring indicates 
the expected dust sublimation radius of 0.05 AU predicted by the size-luminosity-relation (see Sect. 4). 
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ally from 0. 1 1 AU to 0. 14 AU and has a temperature of 1690 K 
at the inner ring edge (see parameters in Table 4). 

4. Discussion 

For comparison with other pre-main sequence stars, we plot the 
geometric /T-band ring-fit radii, where r r i ng j n » 0.13 AU for the 
geometric model without a halo and ^0. 1 1 AU for the model 
with a halo; Sect. 3.1, of S CrA N in the size-luminosity dia- 
gram that shows the /f-band ring-fit radius as a function of the 
stellar luminosity L* (Monnier & Millan-Gabet 2002, Fig. 5). 
The figure shows a sample of Herbig Ae stars (filled squares, 
Monnier et al. 2005), for which this correlation has been orig- 
inally found. Additionally, we plot a sample of TTS (open 
squares, Pinte et al. 2008). 

Figure 5 shows that our measurements of the /f-band ring- 
fit radius of S CrA N of ^0.11-0.13 AU, radii derived with 
the geometric ring-star and ring-star-halo models, is approx- 
imately 2.4 times larger than the dust sublimation radius of 
~0.05 AU predicted for the silicate dust sublimation tempera- 
ture of 1500 K and gray dust opacities (Monnier & Millan-Gabet 
2002). However, several effects can influence the inner model 
radius, such as the chemistry and grain size of the dust (e.g. 
Monnier & Millan-Gabet 2002) or magnetospherical disk trun- 
cation (Eisner et al. 2007). Schegerer et al. (2009) modeled their 
mid-infrared MIDI data and the SED of S CrA N using the 
Monte Carlo code MC3D (Wolfetal. 1999; Schegerer et al. 
2008). They found that an assumed sublimation radius of 
0.05 AU agrees with their observations. 

Furthermore, we compared our ring-fit radii of «0.11- 
0.13 AU derived from the geometric star-disk and star-disk- 
halo models and the inner disk radii of our temperature- 
gradient models (0. 1 1-0. 14 AU) with the predictions of a model 
that accounts for both backwarming and accretion luminosity 




U/L s 

Fig. 5. Location of S CrA N in the size-luminosity diagram: 
K-band ring-fit radii derived with a geometric ring-star model 
(pink square, ^0.13 AU) and a ring-star-halo model (blue square, 
x:0.11 AU); see Tables 3 and 1 for the stellar parameters used. 
The lines indicate the predicted theoretical dependence of the K- 
band radius on the luminosity (Monnier & Millan-Gabet 2002). 
Predictions of models including backwarming are discussed in 
Sect. 4. The filled black squares are a sample of Herbig Ae stars 
(Monnier et al. 2005). The open black squares represent a sam- 
ple of TTS (Pinte et al. 2008). 



(Millan-Gabet et al. 2007; Dullemond & Monnier 2010). This 
model suggests an inner disk radius of 0. 12 AU for a stellar lumi- 
nosity of 2.3 L Q and an accretion luminosity of 0.6 L (derived 
from Prato et al. 2003; Muzerolle et al. 1998), which agrees with 
our measured inner radius of »0. 11-0.13 AU. 

5. Conclusions 

We have observed the T Tauri star S CrA N with AMBER in 
the H and K bands. We have fitted the geometric star-disk and 
star-disk-halo models to our visibility data and derived a radius 
fring.in of approximately 0.11-0.13 AU in the K band (0.095- 
0.12 AU in the //band). We compared the position of S CrA N 
in the size-luminosity diagram with the position of other YSO, 
and found that it is above the line expected for a dust sublima- 
tion temperature of 1500 K and gray dust, but within the region 
of other TTS. The radius predicted by this size-luminosity rela- 
tion is «0.05 AU, whereas the derived ring-fit radius is «0.11- 
0.13 AU (Table 3). However, models including backwarming 
(Millan-Gabet et al. 2007; Dullemond & Monnier 2010) suggest 
a larger inner disk radius of «0.12 AU, which agrees with our 
derived ring-fit radiii of w0. 1 1-0. 13 AU. 

We tested several temperature-gradient models (one- and 
two-component disk models, with or without halo). We found 
that the near- and mid-IR visibilities, as well as the SED, can 
approximately be reproduced by a temperature-gradient model 
consisting of a two-component ring-shaped disk and an un- 
resolved star. The favored temperature-gradient model C (see 
Table 4) has a temperature of »1700 K at the inner disk radius 
of 0.11 AU. The temperature power-law index qi of the inner 
narrow temperature-gradient disk is approximately 0.2, and the 
index qi of the extended outer disk is approximately 0.5, which 
suggests a flared irradiated disk structure (Chiang & Goldreich 
1997). However, q\ is not well-constrained because of the nar- 
row width of the inner component of only s;0.03 AU. The in- 
ner temperature-gradient disk radius of 0.11 AU is similar to 
the four geometric ring-fit radii of approximately 0.10 AU to 
0.13 AU and to the prediction of models including backwarming 
(0.12 AU). Unfortunately, we cannot place any constraints on 
the gas within the inner disk radius of »0. 1 1 AU since the disk 
is only partially resolved (all NIR visibilities are > 0.77) and 
the visibility errors are large. Interestingly, the inner disk com- 
ponents of all three best-fit temperature-gradient models consist 
of a very narrow, hot ring (T^i ~ 1500-1700 K; ring width only 
^0.03-0.04 AU) surrounded by a colder (<600 K) disk com- 
ponent with an extension of several AU. This size and temper- 
ature structure is similar to the structure of more sophisticated 
radiative transfer models including a hot PUIR (e.g., Natta et al. 
2001 ; Dullemond et al. 2001). This suggests that the derived nar- 
row, inner disk component in the temperature-gradient models is 
caused by a hot, perhaps PUIR, in S CrA N. This rim may ex- 
plain the steep temperature jump from 1500 to 600 K. 
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Appendix A: SED references 

To simplify the fitting process, we did not use the pile of original 
data (Table A. 1 ), but instead binned them to form a smaller num- 
ber of data points. Figure A.l shows both the original data (blue 
crosses) and the binned data (red dots). The measurements were 
not recorded contemporaneously (see Table A.l) and therefore 
variability is one of the error sources. Other SED error sources 
are the different aperture sizes and the small number of observa- 
tions in the visible and NIR. 



Table A.l. Measurements from the literature that were used to 
produce the binned SED of S CrA N. A beam width of 10" or 
larger most probably also includes the light from the secondary 
S CrA S (separation 1.4"). 



Wavelength range 


Reference 


Beam 


2.47-11.62 


ISO spectrum 




3.5-170 


ISO 




53.6-106 


Spitzer MIPS 


2.7 x 0.34' 


5.13-36.9 


Spitzer IRS 


55x81" 


12-100 


IRAS 


10" 


450-1100 


JCMT (upper limit) 


10" 


12-100 


Gezari catalog 


10" 


0.768 


Denis 2005 


0.5" 


0.429, 0.554 


Tycho catalog 


10" 


0.657, 0.429 


USNO 


10" 


0.554 


WDS 


0.5" 


1.25-3.55 


Prato 2003 


0.5" 


2.18-18.1 


McCabe 2005 


0.5" 


1.25-3.55 


Morel 1978 


10" 



io- 1 
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X/fim 



Fig. A.l. SED of S CrA N - original data according to the 
sources in Table A.l (blue crosses) and the binned data (red 
dots). 
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